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Advances in the understanding of the basic properties of QCD vacuum will be 
reported. Three main subjects will be touched: 

1) Condensation of monopoles and confinement. 

2) Topology, or instanton physics. 

3) Gauge invariant field strength correlators, and their behaviour across the 
deconfining phase transition. 

1 Introduction 

I will report on some progress recently achieved mainly by numerical simu- 
lations on the lattice, on three subjects. 1) Condensation of monopoles and 
colour confinement (sect. 2). 2) Topology (instanton physics) (sect. 3). 3) Gauge 
invariant field strength correlators at short distances at T = and at the de- 
confining transition (sect. 4). 

The conclusions of this review will be that: 



1) Solid evidence exists that dual superconductivity Elula is the mechanism 
of colour confinement OB 

2) The Witten-Veneziano ifl formula for the mass of the rf is definitely 
correct. This statement has been made possible by the construction of an 
improved (almost perfect) operator for the topological charge densityp By 
use of the same technical improvement a detailed study of the topological 
susceptibility across the,deconfining transition has been made, showing 
a sharp drop of it at T c a 

Studies with dynamical fermions have been started, especially with the 
aim of determining, in full QCD, Xi its derivative \' with respect to 
momentum transfer, and the matrix element of topological charge on 
proton states, which would allow a measurement from first principles of 
the so called "spin content of the proton" . All these results in full QCD 
have been slowed down by the discovery 113 that the usual algorithm 
for Montecarlo simulation, the so called hybrid montecarlo, is very slow 
in changing the topological charge Q of a configuration: the same Q 
stays for few hundred updatings. This puts severe limitations on the 
validity of usual montecarlo studies in full QCD. A typical statistics is 
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indeed a few hundred configurations. Only for observables which are 
insensitive to topology can such a statistics be sufficient to insure proper 
thermalization. 

3) Previous studies of the field strength correlators in the vacuum have 
been extended to much shorter distances (~ 1 fm) . The results is relevant 
to test stochastic models of confinement. 

All these results have been made possible by the use of QUADRIX computers. 



2 Condensation of monopoles and confinement. 

An appealing mechanism for colour confinement in QCD is dual supercon- 
ductivity of the vacuum. The chromoelectric field acting between heavy q q 
pairs is constrained into Abrikosov flux tubes, generating a potential energy 
proportional to the distance. 

We teat this mechanism from first principles by numerical simulation on 
the latticecm 

The basic idea of the approach is that, to detect dual superconductivity of 
QCD vacuum, the vacuum expectation value of an operator with non trivial 
magnetic charge can be used as a probe, or as disorder parameter. Indeed a non 
zero value of such v.e.v. signals condensation of monopoles and spontaneous 
breaking of the U(l) symmetry related to magnetic charge conservation. The 
relevant magnetic f7(l)'s are identified by the well known procedure of abelian 
projectionl!3 They correspond to residual U(l) invariance after diagonalization 
of any operator belonging to the adjoint representation of the gauge group. 
Each operator in the adjoint representation defines an abelian projection, and 
hence a species of monopoles. 

We have constructed a creation operator for monopole, and we use its 
v.e.v. as a disorder parameter. 

It is an open question if all abelian projections are equivalent, and identify 
dual superconductorsEj We have tested the abelian projection defined by the 
diagonalization of the Polyakov line, which is the local operator defined as 
the parallel transport on the closed path along the time axis, starting from a 
point and coming back to it through the periodic boundary conditions used 
to define finite temperature. The v.e.v. of the trace of such operator is the 
usual order parameter for confinement. We find definite evidence for dual 
superconductivity in the U(l) defined by such abelian projection^ 

The creation operator for an U(l) monopole is defined as follows a: 



m(£ t) = cxp 



d 3 xE(x,t)-b(x- y) 
e 



(1) 
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E is the electric field operator, and 

b(f) = ■ 



2 r(r — r ■ n) 

is the vector potential describing the field of a monopole of magnetic charge m 
in units l/2e. The prescription for the integral is that the Dirac string must 
be removed, p is the analog of the translation operator for a particle: 

e ipa \x) = \x + a) 

E is the conjugate momentum to the field A. In the Schrodinger representation 

Li(y,t)\A(x,t)) = \A(2,t) + -b{x - y)) (2) 

e 

/i adds a monopole to the field configuration. Moreover, if Qm = J d 3 x V(V A 
A) is the magnetic charge 

[Qm ,fj,(y,t)]=m(i{y,t) (3) 

i.e. \x carries magnetic charge m. 
The lattice version of ll is 

H = cxp - ((3S M ) (4) 

with 

5M = ^n 0l (n,n )(e it --l) (5) 

ft 

ilmlloi is the lattice version of the electric field. The disorder parameter is 



dUexp[-/3(S w +S M )\ 

(M> = J ~ t (6) 

/ dUexp{~f3S w ) 

Eq. H corresponds to Eq. |l| at the lowest order in b, but it is a compactified 
version of it. By use of the definition Eq. the gauge arbitrariness in the choice 
of b to describe the monopole field is reabsorbed in the Haar measure of the 
Feynman integral. To avoid problems with fluctuations we do not measure (/z) 
directly, but p — 4^ In ii and reconstruct (/i) as 



(fi) = exp / p{x) dx (7) 



o 



(p) determined in this way converges to a nonzero constant as V — > oo for 
T < T c and goes to zero as (T c — T) s at the deconfining temperature. The 
typical form of p is shown in fig. 1 for 51/(3) on a 12 3 x 4 lattice. The sharp 
negative peak at T c signals the drop of (p) at T c . (p) is..z£ro for T > T c only 
in the limit V — > oo, as is typical of disorder parameterstaHj 
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Fig.l SU(3) gauge theory. Lattice 12 3 x 4. p = ^ ln(p) 

We have repeatedQ the construction for the XY model in 3-d, both as a 
test of the method and as an amusing application of it. The model describes 
superfluid He^. It can be viewed as the euclidean version of a free massless 
particle in 2 + 1 dimension. 

The field variable is an angle 9(x). On the lattice: 



5 = /3^[l-cos(A^)] oSo |(A 



(8) 



The theory admits soliton configurations which have the topology of vortices 



9{x — y) = arctg- 



(9) 



'{x- y)i 

is an example. Vortices condense in the ground state in the superfluid phase: 
they play a similar role as monopoles in QCD. As for QCD we can define a 
creation operator for a vortex 



Mz7; *) = ex P 



d 6(x, t)0(x - y) & 2 x 



(10) 



4 



The operator which counts vortices, V can be written as 



V = — I d 2 x 
2tt 



A = ve 



(ii) 



and 



[V,l4=H (12) 

On the lattice by the same procedure of compactification used for monopoles 



v(y,Xo) = exp { -j3 



^2 cos [ A oO {ft} + 6(n - y)] - cos [A 6(n)] 

n=£0 



(13) 



As in the case of monopoles we define (/i) in terms of the correlator of a vortex 
antivortex by cluster property 



lim (fi(xM0)) ~ Ae- a W + ( M > 2 

\x\ — *-oo 



and instead of (//) we measure 



(14) 



(15) 



The behaviour of p is shown in Fig. 2. 

A finite size scaling analysis, described in refta gives, for the critical index 
of the correlation length 



^ = 0.669 ± 0.065 (0.670(7)) 
for the critical temperature 

p c = 0.4538 ± 0.0003 (0.45419(2)) 
and for the critical index of (/x) 

6 = 0.740 ±0.029 



(16) 



(17) 



(18) 



The first two quantities agree with the determinations by other methods, which 
are shown in parentheses, showing that the method is correct and effective. 
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Fig.2 XY model. Lattice 40 3 . p = Aln(/i) v.s. /3. 



3 Topology (instantons) 

Almost all of the Gell-Mann current algebra, which was originally abstracted 
from massless free quark model, has been preserved with the advent of QCD, 
except for the singlet axial current, j'g = ipfj^^ipf, which was conserved 
in the quark model, but is anomalous in QCD 

« = 2N f Q (19) 

Q( x ) = S^ etll,paG lvG a pa is the density of topological chargeElQ = / d 4 xQ(x) 
is the second Chern number, is known as topological charge and takes on in- 
teger values on smooth configurations. The non consevation of in QCD 
gives a handle to solve the so called U(l) problem: if were conserved and 
the corresponding symmetry were a Wigner symmetry, parity doublets should 
exist, if it were realized a la Goldstone then one siupuld have < ^fim^ and 
neither of these possibilities is realized in natureES 

However in an expansion in 1/N C Q{x) is non leading, being oc g 2 = \/N c 
(A = g 2 N c ). In a philosophy m r 2zhich the leading order in 1/N C describes 
the essentials of hadron physics(3E3 U(l) can be considered a simmetry. The 
anomaly acts there as a perturbation and displaces the rj' mass from .zero, 
which would correspond to Goldstone symmetry, to its true valueEJi23 The 



0.0 
-500.0 
-1000.0 



G 



quantitative relation isE9 

?-J-X = ™% + ml, -2m 2 K orx= (180 MeV) 4 (20) 

Jit 

X is the topological susceptibility of the vacuum 

X = J d 4 x(0\T(Q(x)Q(Q))\0) (21) 
at the leading order, i.e. in the absence of fermions. 



The prediction Eq. 20 can be tested on lattice and involves the computation 
of x m the quenched approximation. The correct way to compute x is to 
define a regularized version of Q, Ql, and compute the lattice topological 
susceptibility 

XL = ^(Qi(n)Qi(O)) (22) 

n 

Like in any regularization scheme the regularized operator mixes, in the limit 
in which the cutoff is removed, with all the operators having the same quantum 
numbers and smaller or equal dimension in masscja 
By use of this prescription one gets 

Ql = Z([3)Q (23) 
XL = Z 2 (f3) X a 4 + B(P)G 2 a i + P(p) (24) 

In Eq. |J G2 = (2* G^GfiJ) is the gluon condensate, and P((3) describes the 
mixing to the identity operator, which is usually called perturbative tail. Both 
the last two terms in Eq. [24| come from the singularity in the definition (^2|) of 
X as x — > 0. Z{(3), B((3), P{(3) depend on the choice of Ql, which is arbitrary 
for terms of 0(a 6 ) or higher. For the most simple choice of Ql, in terms of 
the plaquette IF" 



I 

32^ 2 



^ = WTr {IP»IF»} (25) 



Z = 0.18 and the correction to continuum in Eq. g4|, namely the sum of the 
last two terms in Eq. |^| is much bigger than the first term. 

A nonperturbative determination both of 7, a nd of the additive renormal- 
ization M(j3) = B{(3)G 2 a i + P(f3) is possibleBH The idea is based on the fact 
that changing the number of instantons by local updating procedure is much 
slower process than thermalizing the short range fluctuations which are respon- 
sible for renormalizations. So, starting from a zero field configuration, where 
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Q = and x = 0, and heating it will only produce M(/3). Similarly putting 
one instanton of charge Q on a lattice and heating it will leave Q unchanged 
for a large number of sweeps, but will produce the fluctuations necessary to 
build up Z: a plateau will be reached where Ql = ZQ, and Z can be read on 
it. 

In any case it is unpleasant that most of the observed signal is artifact to 
be removed. Recently, playing on the arbitrariness by higher terms in a in the 
definition of Ql } an improved operator has been constructed for which lattice 
artifacts are reduced by 2 orders of magnitude!! Z is now of order 0.6, instead 
of 0.18, and M((3) is less than 10% of the whole signals in a wide scaling 
window. 

Counterterms are removed by the same non perturbative technique, but 
now they are a small part of xl- In particular P(j3) is negligible in a wide 
range of /3's (~ 2% on the entire signal). 
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Fig. 3 Determination of x with the naive operator, Eq. (0 smear) 
and with the once and twice improved operators. 

The results are summarized in Fig. 3 and in Fig. 4. Fig. 3 shows the 
determination of x a t T — 0: the result is x = 175 ± 5 MeV in excellent agree- 
ment with the Witten Veneziano prediction, and in agreement with previous 
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determinations^! within their large errors. A new result is the behaviour of x 
across the deconfining transition shown in Fig. 5, which is .certainly of interest 
for the models of QCD vacuum based on instanton liquidcj 
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Fig. 4 M (J3) j a 4 (/?) . The slight deviations from a constant allow to extract 
P{(3). The result is P(f3) < 2-3% X l 
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Fig. 5 x across T c 



9 



4 Gauge invariant field correlators. 



The correlators are defined as 

ZW0»0 = (°I T (G^(x)S(x)G pa (0)S\x)) |0) (26) 

with G^v — '^Z a T a G ( l lu {x) (T a generators of the gauge group) and S(x) the 
parallel transport from to a; 

S(x) = cxp ^i J A^(tx)x^d?j (27) 

A, = EaT a A^(x). 

These correlators are of interest for stochastic models of confinement, 
where ttipy play a fundamental role, being the lowest order in a cluster ex- 
pansion EOca of the correlation functions. 

The correlators had been determined on lattice a few years agol^ in the 
range of distances from 0.4 to 1 fin. 

Measuring T> pupa (x) on the lattice naively, taking for G pv the open pla- 
quette oerator, is difficult because of lattice artefacts due to short range fluc- 
tuations. 

If these fluctuations are smeared, e.g. by a local cooling procedure, a 
plateau will eventually be reached in the cooling process, where short range 
effects have been removed, but large distance physics is left unchanged. The 
basic idea is that in a local cooling precess a distance d is reached after a 
number of steps t which is governed by a sort of diffusion equation 

t cx d 2 



The minimum distance which can be explored by this technique is a few 
(~ 4) lattice spacings. 

In physical units the limitation comes from the requirement that the lattice 
size La must be larger than the typical scale of 1 fm. At (3 ~ 6.1 1 fm is 8 
lattice spacings so that distances from 0.5 to 1 fin can be explored on a 16 4 
lattice. This was the range of distances in refO. To have a distance of 3-4 
lattice spacing corresponding to 0.2 fm one needs (3 ~ 7 and a lattice size 
L > 32. FkrJl shows the two independent form factors D and D\ which 
parametrize EZIE3 T> pvpa as a function of distance, extracted from the joint 
sample of data of refill and 29. 

Across T c the magnetic correlators are the same, as at T — 0, while the 
electric correlators drop by an order of magnitudeESa This is consistent with 
the vanishing of the string tension at T > T c . 
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(lower curve). 
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